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Abstract. Following up on earlier work on this topic ( Reinecke et al. 1999b a|), we present an improved set of 
numerical models for simulations of white dwarfs exploding as Type la supernovae (SNe la). Two-dimensional 
simulations were used to test the reliability and numerical robustness of these algorithms; the results indicate that 
integral quantities like the total energy release are insensitive to changes of the grid resolution (above a certain 
threshold), which was not the case for our former code. 

The models were further enhanced to allow fully three-dimensional simulations of SNe la. A direct comparison of 
a 2D and a 3D calculation with identical initial conditions shows that the explosion is considerably more energetic 
in three dimensions; this is most likely caused by the assumption of axisymmetry in 2D, which inhibits the growth 
of flame instabilities in the azimuthal direction and thereby decreases the flame surface. 

Key words, supernovae: general - physical data and processes: hydrodynamics - turbulence - nuclear reactions, 
nucleosynthesis, abundances - methods: numerical 



1. Introduction 

The project described in this ar ticle is a direct continu - 
ation of the work presented by Reinecke et al. (1999b| ), 
and is therefore concerned with the accurate modeling of 
the rapid thermonuclear combustion processes which oc- 
cur in a white dwarf during its explosion as a Type la 
supernova. While the earlier results already gave several 
new insights into the explosion dynamics of such an event 
and its dependence on parameters like the ignition con- 
ditions, it nevertheless became clear that the numerical 
models employed at that time were not sufficient to re- 
produce characteristic properties of SNe la that are known 
from observations; most importantly, the total energy re- 
lease predicted by the simulations lay consistently below 
the expectations. Meanwhile, several possible reasons for 
this deficiency could be identified and the models were 
enhanced and adjusted accordingly, leading to much more 
realistic results. 

The following section gives details of the implemented 
changes and explains why they were required for the cor- 
rect treatment of SNe la. These modifications affect the 
numerical description of the flame as well as the model 
for the unresolved turbulent motions and the thermonu- 
clear reaction rates. Section || presents the results obtained 
from a two-dimensional resolution study, which is mainly 
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intended to test the correctness and numerical robustness 
of the code, even though data obtained from 2D simu- 
lations is unlikely to agree with observations for reasons 
given in section 

The first three-dimensional simulation performed with 
the new code is discussed in section fi], with a focus on 
analyzing the similarities and differences compared to the 
axisymmetric 2D model. 

Finally, in section || our data are compared with obser- 
vations and numerical simulations performed by other au- 
thors; furthermore, several possible future improvements 
of the code are presented briefly. 

Forthcoming articles will focus on a more detailed anal- 
ysis of three-dimensional simulations with varying initial 
conditions and numerical resolution. 



2. Improvements of the numerical models 

It must be kept in mind that the corrections and addi- 
tional features introduced below aim at a more realistic 
description of the explosion dynamics and energetics. No 
effort was made to reproduce other aspects, like the exact 
chemical composition of the burned material, unless this 
was required for the goals mentioned above. As a conse- 
quence, the treatment of nuclear reactions might appear 
minimalistic in comparison to the large reaction networks 
typically used in one-dimensional simulations; but reach- 
ing a comparable level of detail in multidimensions would 
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be prohibitively expensive and is fortunately not required 
for reliable predictions of the explosion strength. 

2.1. Nuclear reactions 

In a SN la, the nuclear reactions can be distinguished 
depending on their time scales: 

— The "fast" reactions, which transform the initial C/O 
mixture into nuclear statistic equilibium (NSE) at high 
densities or into a mixture of intermediate-mass ele- 
ments at lower densities, take place at time scales Tf ast 
that are much shorter than the time T cross during which 
the flame crosses a grid cell: 



Tfast < Tc: 



A/s, 



(1) 



where A denotes the linear extent of a cell and s the 
flame propagation speed. Therefore the transition be- 
tween unburned and burned material is too abrupt to 
be spatially resolved and should be modeled as a nu- 
merical discontinuity. This can be done by using the 
so-called level set technique (see section 2.2). 

— Composition changes due to the expansion of the star 
(i.e. changes of the NSE), on the other hand, are much 
slower than r cross and the Courant time step and can 
be modeled by simply changing the local species con- 
centrations to the desired values after each time step. 

This picture is not entirely correct, since there exist - 
especially at lower densities - reactions behind the flame 
whose time scales are comparable to T cross . However, these 
reactions have no significant effect on the internal energy 
or the molecular weight of the material and therefore can 
be neglected if one is only interested in the hydrodynamics 
of the explosion. 

The original version of the combustion algorithm (see 
Reinecke et al. 1999b) did not contain any slow reactions, 
and the fast reactions were approximated rather crudely 
by the instant fusion of carbon and oxygen to nickel. In 
the meantime we realized that such a description is over- 
simplified and cannot lead to accurate results in a Type 
la model for several reasons: 

— Combustion at high densities does not immediately 
produce nickel, but also a considerable fraction of a- 
particles. The specific energy release of this reaction is 
much lower than for the production of pure nickel. 

— However, as the density in the ashes drops due to the 
stellar bulk expansion, the equilibrium shifts towards 
the heavy elements and the energy "stored" in the en- 
particles is released during the later explosion stages. 

— As soon as the flame enters regions of intermedi- 
ate density (10 8 -10 7 g/cm 3 ), elements of intermediate 
mass are synthesized instead of nickel, and the specific 
energy release drops accordingly. 

Since all of these formerly neglected phenomena may 
affect the explosion dynamics, they were implemented in 
our new scheme, which is described below: 



The initial mixture consists of 12 C and 16 at low tem- 
peratures. Because of the electron degeneracy the fuel 
temperature is nearly decoupled from the rest of the 
thermodynamical quantities, and since temperature is 
not used to determine the initial reaction rates, its ex- 
act value is unimportant. 

When the flame passes through the fuel, carbon and 
oxygen are converted to ash, which has different com- 
positions depending on the density of the unburned 
material. At high densities (p > 5.25 • 10 7 g/cm 3 ), a 
mixture of 56 Ni and a-particles in nuclear statistic 
equilibrium is synthesized. Below that density burn- 
ing only produces intermediate mass elements, which 
are represented by 24 Mg. Once the density drops below 
10 7 g/cm 3 , no burning takes place. 
Because NSE is assumed in the ashes, the proportion of 
56 Ni and a-particles will change depending on density 
and temperature. To obtain the correct values, a NSE 
data table provided by H.-Th. Janka and an iterative 
root finding algorithm were employed. 
Whenever the composition in a cell changes due to 
nuclear reactions, a corresponding amount of energy 
is added or subtracted from the total energy in that 
cell. The amount is d erived from the n uclear binding 
energies published by Audi et al. (1997| ). 



Compared to the much simpler original approach, this 
more detailed and realistic model has the potential to al- 
ter the explosion dynamics significantly: first, the reaction 
to NSE and intermediate elements with its comparatively 
low heat release results in lower temperature and higher 
density of the burned material, which will decrease the 
buoyancy of the burning bubbles. In addition the number 
of particles per unit mass in the ashes is also higher with 
the new model, which reduces the temperature even more. 
Since the a-particles will be converted to nickel during the 
expansion of the star, not all of the nuclear energy is re- 
leased at once when material is processed by the flame, 
but the release is partially delayed. The combination of 
all these effects is expected to shift the maximum of en- 
ergy generation to a later time during the explosion. 

The transition densities from burning to NSE and in- 
complete burning, as well as from incomplete burning to 
flame extinction were derived from data of a W7 run pro- 
vided by K. Nomoto. This approach is rather phenomeno- 
logical, and since these densities can have a potentially 
large impact on the simulation outcome, it will have to be 
re-examined in a thorough manner. 



2.2. Thin flame model 

The numerical representation of the thermonuclear reac- 
tion front (i.e. the location where the "fast" reactions take 
place) did not change f undamentally compare d to the orig- 
inal implementation ( Reinecke et al. 1999b ). To repeat 
briefly, the flame front is associated with the zero level set 
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of a function G(r, t), whose temporal evolution is given by 



dG 
dt 



= -(v u + s u n)(-n\VG\ 



(2) 



where v u and s u denote the fluid and flame propagation 
velocity in the unburned material ahead of the front, and 
n is the front normal pointing towards the fuel. The ad- 
vection of G caused by the fluid motions is treated by the 
piecewise parabolic method ( |Colella fc Woodward 1984 ), 
which is also used by our code to integrate the Euler equa- 
tions. After each time step, the front is additionally ad- 
vanced by s u At normal to itself. 

This equation is only applied in the close vicinity of 
the front, whereas in the other regions G is adjusted such 
that 



IVGI = 1. 



(3) 



The source terms for energy and composition due to 
the fast thermonuclear reactions in every grid cell are de- 
termin pd a3 follow3 i 



In our simulations, l g almost always lies far below the grid 
scale. Simply neglecting the surface increase on sub-grid 
scales would lead to an underestimation of the energy gen- 
eration rate, which is not acceptable; therefore a model for 
a turbulent flame speed St > s/ is required to compensate 
this effect. 

For the case of very strong turbulence (i.e. for v'(A) ^> 
si or l g <C A for a numerical resolution of A) it has 
been shown that the turbulent flame velocity decouples 
fro m si and is pro p ortional to the velocity fluc tuations 
v' dShchelkin 19431 ; [Pocheau 19921 ; [Peters 2000|) . In our 
simulations this condition is fulfilled except at the very 
beginning of the explosion, where combustion starts lam- 
inar before sufficient turbulence has been generated; this 
stage only lasts a few tens of milliseconds. At lower den- 
sities (p J$ 10 7 g/cm 3 ) this picture is possibly not correct, 
since in that case the combustion takes place in the so- 
called distributed burning regime (Peters 1986); it is still 
a matter of debate whether a transition to a detonation 



can occur at these densities (Niemeyer 199S; Khokhlov 



et al. 1997), which would of course lead to a quite dif- 



^Ashes = max(l - a, A Ashes ) 
= 1 



^Fuel 



'tot 



etot 



Y' 

Ashes 



Ashes 



x 



Ashes 



(4) 
(5) 
(6) 



where a is the volume fraction of the cell occupied by un- 
burned material; this quantity can be determined from the 
values of G in the cell and its neighbours. The quantity q 
represents the specific energy release of the total reaction. 

In the new model for the nuclear reactions (see section 
2.1) the composition of the ashes depends on the density 
of the unburned material ahead of the front, which cannot 
be properly determined from the state variables of the grid 
cell in question, since it contains a mixture of burned and 
unburned states from which a reconstruction is not possi- 
ble. Currently the maximum density of all neighbour cells 
is adopted as an approximation for the real density of the 
unburned material and used to determine the compositon 
of the newly created ashes. This approach is acceptable 
in the context of SNe la, since the white dwarf does not 
contain steep density gradients (with the exception of the 
flame itself) , and supersonic phenomena, which might cre- 
ate such gradients, do not occur. 

2.3. Model for the turbulent flame speed 

All multidimensional simulations of exploding white 
dwarfs share the problem that it is impossible to resolve 
all hydrodynamically unstable scales. The consequence is 
that the simulated thermonuclear flame can only develop 
structures on the resolved macroscopic scales, while the 
real reaction front will be folded and wrinkled on scales 
down to the Gibson scale l g , which is indirectly defined by 



ferent chemical composition and a higher energy release 
in the outer stellar layers. For the time being we ignore 
the possibility of a delayed detonation; if further work in 
this field indicates that this scenario cannot be ruled out, 
however, this point will have to be addressed. 

As discussed above, the quantity v'(A) must be known 
in order to determine the turbulent macroscopic flame 
speed St- Unfortunately it cannot be determined directly 
from the Reynolds stress tensor at the location of the 
front because of the velocity jump between burned and 
unburned material. Therefore we make use of a technique 



first presented by Clement (1993) and later applied to SNe 
la by Niemeyer & Hillebrandt (199£), which models the 
creation, advection and dissipation of the turbulent kinetic 
energy on sub-grid scales q from which v'(A) can be easily 
derived. 

The time evolution of q is given by 



d{pq) , 2 __ 

+ V{vpq) = --pqVv 



V(^turbVg). 

(8) 



The individual source and sink terms on the right hand 
side can be interpreted as turbulent compression, input 
from macroscopic scales, dissipation into thermal energy 
and turbulent diffusion, respectively. 

Following Clement (1993 ), the turbulent viscosity, 



stress tensor and energy dissipation are modeled by 



rb = pAv'(A) = pCA^, 



v'(l g ) = si, 



(7) 



and 



where v'(l) denotes the amplitude of velocity fluctuations 
on a scale I and s/ the laminar flame propagation speed. 



^dis 



t/ 3 A~ 



= Vq a l 2 A-\ 



(9) 



(10) 



(11) 
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The remaining free parameters C and V are set to 

C = O.LT 7 and (12) 
V = 0.5J 7 -\ (13) 
where T = min(100, max(0.1, lCT^/g)), (14) 

again in accordance to Clement's work. 

From the appearance of only one cell dimension A in 
the above equations it is obvious that this approach is 
only valid for grids with quadratic and cubic cells. For 
anisotropic grids, Clement suggests to calculate ryt ur b sep- 
arately for each spatial direction, but to treat the dissipa- 
tion still as an isotropic effect. The expressions for rfturb 
and ejisg then become 



?7turb,j = pCAf A 1 ^fq and 
e diss =2V /2 A- 1 , 



(15) 
(16) 



where A is the geometric mean of the individual cell di- 
mensions. However, even if this extension is used, the 
model is expected to perform rather poorly in the ex- 
treme case of very elongated cells, since the concept of 
a "cell size", on which q is defined, breaks down under 
these circumstances. 

The model described above produces an approxima- 
tion of the turbulent sub-grid energy throughout the com- 
putational domain, although the source and sink terms for 
q must be disabled in the cells cut by the front because of 
the velocity jump. This is possible because q is transported 
into these cells by means of advection and turbulent dif- 
fusion. The desired magnitude of the velocity fluctuations 
is then given by 



./(A) = y&z, 



(17) 



which completes the model for the flame propagation 
speed. 

In the given form, the sub-grid model is correct for 
three-dimensional simulations only; in two spatial dimen- 
sions it is necessary to replace the factor 2/3 in equa- 
tions (||) and (|Io| ) by a factor 1 to preserve the traceless- 
ness of Sy. Leaving this factor unchanged would imply 
three-dimensional turbulent motions on unresolved scales, 
whereas the flow on the macroscopic scales is only two- 
dimensional. Since the transition scale between both re- 
gions is identical to the resolution, one would expect dif- 
ferent simulation results for identical initial conditions, if 
the cell sizes are varied; this is certainly not desirable. 
In the older versions of the simulation code this fact was 
overlooked, and all two-dimensional calculations were per- 
formed with the incorrect factor 2/3; this resulted in a 
rather strong dependence of the total energy release from 
the numerical resolution. With the corrected model, how- 
ever, the resolution does not influence the explosion ener- 
getics significantly (see section . 

It must be emphasized at this point that any attempt 
to simulate turbulent flow in two dimensions will likely 
produce results that differ significantly from those of fully 
three-dimensional simulations; this effect is caused by the 



unrealistic absence of turbulent velocity fluctuations along 
the third coordinate axis, which may have a strong influ- 
ence on the scaling behaviour of the turbulent motions. 
Therefore the results of two-dimensional calculations in- 
volving turbulent flow must be checked at least punctually 
against their three-dimensional equivalents in order to give 
them the required credibility. 

3. Two-dimensional calculations 

The goal of the calculations presented here was to ver- 
ify that variations of parameters without direct physical 
relevance (e.g. numerical resolution) do not affect the re- 
sults significantly. Because of the enormous computing 
and storage requirements for this kind of numerical study, 
the simulations had to be performed in two dimensions 
assuming axisymmetry instead of solving the full three- 
dimensional problem. Nevertheless, the insights gained 
from these experiments can - at least in a qualitative man- 
ner - also be applied to the three-dimensional case- 



in analog y to the calculations described in [Reinecke 
|et al. (1999a ), the experiments were carried out in cylin- 
drical (r, z) coordinates with a uniform grid spacing in 
the inner regions of the computational domain (from the 
origin to a radius of » 2.3- 10 8 cm). In the outer regions the 
radial and axial dimensions of the grid cells were increased 
exponentially in order to avoid mass loss across the grid 
borders due to the stellar expansion. Equatorial symme- 
try was assumed; at the same time rotational symmetry 
around the polar axis was imposed by the choice of coor- 
dinates. The white dwarf was constructed by assuming a 
central density of 2.9- 10 9 g/cm 3 , a uniform temperature of 
5 • 10 5 K and a composition of Xq = Xo = 0.5, resulting 
in a mass of 2.797 • 10 33 g, a radius of about 1.8 • 10 8 cm 
and a binding energy of 5.19 • 10 50 erg. 



3.1. Resolution study 

To study the robustness of our code with respect to a 
change of the numerical resolution, simulations were per- 
formed with grid sizes of 128 2 , 256 2 , 512 2 and 1024 2 cells, 
whose corresponding resolutions in the uniform inner part 
of the grid were 2-10 6 cm, 10 6 cm, 5-10 5 cm and 2.5-10 5 cm. 
The initial flame geometry (called c3_2d) used for all 
the se calculations is ident ical to the setup C3 presented 
by Reinecke et al. (1999a ): the matter within a radius of 
1.5 • 10 7 cm from the stellar center was incinerated, and the 
surface of the burned region was perturbed to accelerate 
the development of Rayleigh- Taylor instabilities. 

Figure |l| shows the energy release of the models; ex- 
cept for the run with the lowest resolution, the curves are 
nearly identical in the early and intermediate explosion 
stages. Simulation c3_2d_128 exhibits a very slow initial 
energy increase and does not reach the same final level 
as the other models. Most likely this is due to insufficient 
resolution, which leads to a very coarsely discretized ini- 
tial front geometry and thereby to an underestimation of 
the flame surface. From this result it can be deduced that 
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Fig. 2. Comparison of two simulations with identical initial conditions, but different resolution. The front geometry is 
indicated by a black line, whereas the amplitude of turbulent velocity fluctuations is given color-coded from (white) 
to 2.5T0 7 cm/s (black). The low-resolution run clearly exhibits less flame structure, which is compensated by a higher 
flame speed. The overall features appear to be remarkably similar in both diagrams. 



Resolution study: energy release 
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Fig. 1. Time evolution of the total energy for the initial 
flame geometry c3_2d and different resolutions. During 
the early and intermediate explosion stages there is excel- 
lent agreement between the better resolved simulations. 



all supernova simulations performed with our code should 
have a central resolution of 10 6 cm or better. 

In the late explosion phase (after about 0.8 seconds) 
the total energy is no longer equal for the three simu- 
lations, but increases with better resolution. The origin 
of this effect is not clear, but since it occurs at the time 
when the flame enters the anisotropic part of the grid, it is 
probably caused by inaccuracies of the burning and hydro- 
dynamic algorithms in the highly elongated cells in these 
regions. A breakdown of the sub-grid turbulence model is 
the most likely culprit, since in a very anisotropic grid the 
concept of a "grid scale" is not well defined anymore. In 
any case the observed scatter in the final energy releases 
of the order of 10% was assumed to be acceptable. 

Overall, our model for the turbulent flame speed ap- 
pears to compensate the lack of small structures in the 
front very well. 



Another argument for the reliability of the turbulence 
model is provided by figure [|, which shows the front ge- 
ometry and turbulence intensity after one second for the 
simulations c3_2d_256 and c3_2d_1024. These data are 
in very good agreement with the expectations: 

- Both simulations have quite similar large-scale fea- 
tures; e.g. the burned areas appear similar in size, and 
the largest distance from the front to the center also is 
nearly equal. 

- Because of the larger scale range in c3_2d_1024, the 
flame structure is much better resolved: it shows, for 
example, the onset of Kelvin-Helmholtz-mstabilities in 
the shear flow between rising hot material and falling 
fuel near the coordinate axes and also small, secondary 
Rayleigh- Taylor instabilities, which are not visible in 
the coarser simulation. 

- The amplitude of sub-grid fluctuations is lower for the 
better resolved simulation, which is in agreement with 
the turbulent velocity scaling law. This is also required 
to compensate the larger flame surface and keep the 
total energy generation rate resolution-independent. 

This plot also demonstrates convincingly that the high- 
est turbulence intensities are reached in the shear layer 
between fuel and ashes. 



3.2. Comparison to earlier 2D simulations 

Comparison of the energy evolution with the results pre- 
sented by Reinecke et al. (1999a) reveals that the final 
absolute energies are considerably higher for the refined 
code than for the older calculations. The remnant, which 
had formerly remained bound, now reaches a positive total 
energy and will therefore continue its expansion. However, 
the energy release is still too low to account for a typical 
SN la. 

Significant differences also exist in the time evolution 
of the energy generation rate (figure H) : the period of most 
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Energy generation rate Model c3_2d_256: composition evolution 



1.00 



c3_2d_256 

C3 (Reinecke et al. 1999) * 





Fig. 3. Time evolution of the energy generation rate for 
the models C3 in |Reinecke et al. (1999a| ) and c3_2d_265. 
The shift of the maximum is clearly visible. 



intense burning, which was originally reached after 0.2- 
0.5 s, has now shifted to 0.5 -0.8s after ignition. The main 
reason for this delayed energy release is the more realistic 
treatment of the fusion reactions. Since carbon and oxy- 
gen are no longer instantaneously fused to nickel, but to 
a NSE mixture of nickel and a-particles, less energy is 
released in the early burning phases, which results in a 
lower temperature and higher density in the ashes. This 
in turn decreases the buoyancy and rising speed of the 
burning bubbles, leading to less shear and lower turbulent 
burning speeds. As a compensation, the transformation of 
a-particles into nickel injects additional energy into the 
explosion at later times. 

This argument is supported by an evolution diagram of 
the white dwarf's composition, shown in figure [|. As long 
as the nuclear reactions take place near the center and 
bulk expansion is still rather slow, nickel and a-particles 
are produced at a nearly constant ratio, and the specific 
energy release is much lower than for synthesis of pure 
nickel. During the continued expansion, the equilibrium 
concentration for a-particles in the NSE rapidly drops to- 
wards zero and the "buffered" energy is released as thermal 
energy, further driving the expansion. This very exother- 
mic process (Q « 1.52 • 10 18 erg/g) coincides rather well 
with the maximum of the energy generation rate (figure 

!)• 



Fig. 4. Time evolution of the chemical composition for 
model c3 2d 256. 



The initial configuration, as well as snapshots at later 
times, are shown in figure ||. Obviously, the initial axisym- 
metry is lost after 0.2 -0.3 s, although no explicit pertur- 
bation in (^-direction was applied to the front. This hap- 
pens because the initial flame geometry cannot be mapped 
perfectly onto a Cartesian grid and therefore the front is 
not transported at exactly the same speed for all ip. On a 
cylindrical grid, which matches the symmetry of the prob- 
lem setup better, this effect would not be observed. 

In our case, however, this symmetry breaking is de- 
sired, since an exactly axisymmetric calculation would 
only reproduce the results of the 2D run. Besides, perfect 
rotational symmetry would never occur in reality. 

During the next few tenths of a second, the small devi- 
ations cause the formation of fully three-dimensional RT- 
mushrooms, leading to a strong convolution of the flame. 
As expected, this phenomenon has a noticeable influence 
on the explosion energetics; this is illustrated in figure 
[| Before the loss of axial symmetry in c3_3d_256, the 
total energy evolution is almost identical for both simu- 
lations, which strongly suggests that the two- and three- 
dimensional forms of the employed turbulence and level 
set models are consistent, i.e. that no errors were intro- 
duced during the extension of these models to three di- 
mensions. In the later phases the 3D model releases more 
energy as a direct consequence of the surface increase 
shown in figure |H[ 



4. Extension to three dimensions 

In order to compare two- and three-dimensional simula- 
tions directly, a 3D calculation was performed using the 
same initial conditions as given in section [|. For this pur- 
pose the initial two-dimensional flame location was ro- 
tated by 90 degrees around the z-axis and mapped onto 
the three-dimensional Cartesian grid consisting of 256 3 
cells with a central resolution of 10 6 cm. Only one octant 
of the white dwarf was simulated and mirror symmetry 
was assumed with respect to the coordinate planes. 



5. Discussion and conclusions 

Since the explosion energy and the composition of the 
remnant are important criteria for the quality of a SN 
la model, the relevant data for all performed simulations 
are listed in table 0. 

Given the stellar binding energy of 5.19 • 10 50 erg, it is 
evident that the progenitor becomes unbound in all ex- 
periments, which implies that no recontraction will occur 
and no compact object remains. Nevertheless only model 
c3_3d_256 results in a powerful enough explosion to qual- 
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t=0.15s 



model name 


M 


m Ni 

M 


B, lnc 
10&U erg 


c3_2d_128 


0.132 


0.348 


6.56 


c3 2d 256 


0.109 


0.400 


7.19 


c3_2d_512 


0.151 


0.402 


7.58 


c3 2d 1024 


0.152 


0.428 


8.00 


c3 3d 256 


0.177 


0.526 


9.76 



Table 1. Overview over element production and energy 
release of all discussed supernova simulations 



t=0.45s 




t=0.6s 




t=0.79s 




Fig. 5. Snapshots of the flame front for a centrally ignited 
three-dimensional scenario. One ring on the coordinate 
axes corresponds to 10 7 cm. 



2D/3D energy comparison 




1.5 



Fig. 6. Comparison of the explosion energy for identical 
initial conditions and resolution in two and three dimen- 
sions. After the loss of axial symmetry (at t w 0.3 s) the 
larger flame surface in the three-dimensional model leads 
to more vigorous burning. 



ify as a typical SN la; the two-dimensional scenarios are 
too weak to accelerate the ejecta to the speeds observed in 
real events and produce too little nickel to power a stan- 
dard SN la light curve, which had more or less to be ex- 
pected, since the surface growth was partially suppressed. 

The 3D calculation is a good candidate for typical SN 
la explosions, at least with respect to explosion strength 
and remnant composition. Its nickel mass falls well into 
the range of ~ 0.45 - 0.7 Mq determined by Contardo et al, 
(2000| ) for several typical events, and it can be deduced 



from the amount of "magnesium" in the ejecta that enough 
intermediate mass elements were synthesized to explain 
the observed spectral features. 

Qualitatively, our results for the explosion energetics 
are in rath er good agreemen t with recent simulations per- 
formed by Khokhlov (2000| ), which rely on quite differ- 
ent numerical models and assumptions; most notably the 
macroscopic flame velocity is determined by the asymp- 
totic rise speed of a Ray leigh- Taylor bubble with size A, 
instead of the turbulent velocity fluctuations. The par- 
ticular setup described in his paper is centrally ignited, 
and no perturbation is applied to the spherical flame sur- 
face. Due to the high resolution the discretization errors 
are quite small, and there is a rather long phase of very 
slow burning until the first instabilities have reached the 
nonlinear stage; from this moment on, however, the explo- 
sion progresses in a way very similar to the 3D simulation 
discussed here. This similarity applies to the energy pro- 
duction rate as well as the geometrical features developed 
by the flame. 

Further improvements of the SN la simulations are 
planned, both with respect to the underlying models and 
the realism of the initial conditions. As was pointed out 
by S. Blinnikov (personal communication), an accurate 
simulation of light curves and spectra based on the re- 
sults of the 3D hydrodynamical computation requires not 
only the total concentration of iron group elements (which 
have been represented so far by 56 Ni), but the amounts of 
the real 56 Ni and the rest of the Fe-group elements. The 
necessary changes will be implemented in the near future. 

It would also be desirable to reach a higher numerical 
resolution at the beginning of the explosion, which allows 
a more accurate prescription of the flame geometry, and, 
at the same time, to follow the expansion of the star over 
a longer time period (several seconds), until the ballis- 
tic expansion stage has been reached. To achieve both of 
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these goals, an expanding grid must be used instead of the 
currently employed static one. 
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